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Abstract 



The physics of the inflationary universe requires the study of the out of 
equilibrium evolution of quantum fields in curved spacetime. We present the 
evolution for both the geometry and the matter (described by the quantum 
inflaton field) by means of the non-perturbative large N limit combined with 
semi-classical gravitational dynamics including the back-reaction of quantum 
fluctuations self-consistently for a new inflation scenario. We provide a cri- 
terion for the validity of the classical approximation and a full analysis of 
the case in which spinodal quantum fluctuations drive the evolution of the 
scale factor. Under carefully determined conditions, we show that the full 
field equations may be well approximated by those of a single composite field 
which obeys the classical equation of motion in all cases. The de Sitter stage 
is found to be followed by a matter dominated phase. We compute the spec- 
trum of scalar density perturbations and argue that the spinodal instabilities 
are responsible for a 'red' spectrum with more power at longer wavelengths. 
A criterion for the validity of these models is provided and contact with the 
reconstruction program is established. 



I. INTRODUCTION AND MOTIVATION 
A. Preliminaries 

During the last decade Particle Physics Cosmology matured from the realm of specula- 
tion to that of testable predictability. The next generation of satellite and balloon borne 
experiments will provide an unprecedented test of the fundamental ideas of Early Universe 
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Cosmology and will validate or rule out current theories of structure formation and other 
aspects of the cosmology of the early universe Standard Big Bang cosmology is based 

on a homogeneous and isotropic expanding universe and is confirmed by observations: i) 
redshifts of objects far away have been measured and confirm Bubble's law of expansion, 
ii) the cosmic microwave background radiation (CMBR) has been measured to be an al- 
most perfect blackbody with temperature T = 2.728 ± 0.002/^ and is a remnant of the 
early hot and dense stages of evolution after the Big Bang, iii) light element abundance of 
D,^ He,'^ HeJ Li which is consistent with what is predicted by Nucleosynthesis. 

The CMBR gives evidence of homogeneity and isotropy on scales larger than about 100 
Mpc, with small temperature fluctuations AT/T < 10~^ on angular scales that range from 
1° to about 90° This high degree of homogeneity and isotropy on large scales presents 
one of the important puzzles of standard big bang cosmology, i.e. the horizon problem. 

Distances in a spatially fiat, homogeneous and isotropic cosmology are measured with 
the FRW metric 



with t the comoving or cosmic time and a{t) the scale factor. Physical distances along null 
geodesies determine the limit of causal correlations. These are given by 



this distance determines the causal horizon, events that are separated by distances larger 
than dh{t) cannot be causally correlated by microphysical processes. Another important 
physical scale is the Hubble radius dnit) = H~^(t) = {a{t) / a(t))~^ . In the cases of interest, 
dh{t) and dnit) are proportional and we will call dnit) the horizon indistinguishably. Let 
us consider a physical distance today A (to) over which the CMBR determines homogeneity 
and isotropy to one part in 10^. At a time t earlier, the size of this patch is given by 
X(t) = X{to)a{t) (where we have chosen a{to) = 1). For a matter or radiation dominated 
universe (most of the life of the Universe), a{t) oc (n = 1/2 for radiation, n = 2/3 for 
matter domination). The horizon size (or Hubble radius) is given by dnit) oc t and the ratio 
X{t)/dH{t) increases as we evolve back to an earlier epoch. Thus at some time in the far 
past, the scale A(to) will become larger than the horizon. In particular, if \{to) ~ 100 Mpc 
this scale will be larger than the horizon near the time of decoupling of matter and radiation, 
at red-shifts z ~ 1100, and therefore it is extremely difficult to explain how such a region 
which was causally disconnected at decoupling can be so homogeneous and isotropic. This 
is the essence of the horizon problem |1]J^. 

Inflation was proposed over a decade ago to solve long standing problems of Standard Big 
Bang cosmology such as the homogeneity and horizon problems . Inflation corresponds 
to a stage of exponential expansion of the Universe during which the Hubble radius (horizon) 
remains constant and physical scales grow exponentially. Thus, in this scenario, microphys- 
ical scales smaller than the horizon over which causal processes establish correlations cross 
the horizon during the inflationary stage. After the inflationary stage and through particle 
physics processes (reheating), the Universe becomes radiation and matter dominated and 
the scales that had crossed outside the Hubble radius during the inflationary stage, re-enter 



ds'^ = dt^ — o?{t)dx^, 



(1.1) 




(1.2) 
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the horizon. Perturbations on these scales then grow under gravitational (Jeans) instabil- 
ity to form the large scale structures that we see today. Thus the inflationary proposal 
solves the homogeneity and horizon problems by allowing scales to cross the horizon twice, 
from being subhorizon initially becoming superhorizon during inflation, and re-entering the 
horizon during the period of radiation or matter domination. 

Quantum field theory combined with particle physics models provides the framework 
to implement the inflationary idea. For an inflationary stage the equation of state has to 
be dominated by a term similar to a cosmological constant leading to a negative pressure. 
Through the Einstein-Friedmann equations, this yields to an exponential expansion of the 
scale factor. There are many different inflationary scenarios: 'old', 'new', 'natural', 'hybrid' 
etc. [mi , but despite the plethora of different implementations of the original idea, there are 
some robust features of inflation common to practically all models: i) Inflation predicts aflat 
universe, i.e. the total energy density at all times is the critical energy density, ii) an almost 
scale invariant spectrum of scalar density perturbations, with amplitudes that are bound by 
the CMBR inhomogeneities AT/T, iii) approximately scale invariant spectrum of gravita- 
tional waves (tensor perturbations). The next generation of satellites (MAP/PLANCK) will 
provide more stringent bounds on temperature inhomogeneities and a firmer determination 
of the spectrum of scalar and tensor perturbations, and the next generation of gravita- 
tional wave detectors (LIGO, VIRGO, LISA) may detect the gravitational waves and their 
spectrum. Determinations of the spectrum of scalar and tensor perturbations will provide 
stringent bounds on inflationary models based on particle physics scenarios and probably 
will validate or rule out specific proposals. 

Although the general features of inflation seem robust and universal, the implementation 
of the inflationary stage as well as the departure from scale invariance of the spectra of scalar 
and tensor perturbations depend on the concrete model and its dynamics. In this article 
we focus on the study of the dynamics of specific new inflationary scenarios that lead to an 
inflationary stage after a phase transition. 



B. Dynamics of Phase Transitions 

An appealing new inflation scenario envisages a phase transition at GUT scales that 
provides an inflationary stage . The concept is an extension of well understood features 
of phase transitions in model field theories and has been studied within the context of 
inflationary cosmology in a fixed background by Linde and Vilenkin |^,^ and Guth and Pi 
||T0| . Consider a scalar field theory (we will neglect fermions, gauge fields and other exotic 
components in the discussion) with a typical 'Mexican hat' potential that allows for broken 
symmetry states. If the system is originally at very high temperature, larger than critical, 
i.e. Tj > Tc, the 'effective' potential has a minimum when the expectation value of the scalar 
field (the order parameter) vanishes. As the temperature cools below critical, the effective 
potential develops minima away from the origin and the expectation value of the scalar field 
will roll toward the minima, i.e. the equilibrium values. Consider as a relevant example the 
following scalar potential 

A 1 171^ 

Veff{<P, m) = 4 0' + 2 ni'[T{t)] 02 + — , (1.3) 
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m'[T]=m^ S-1 , (1-4) 

T' Tft 

where the time dependence in the temperature results from considering a resummation of 
the fluctuation contribution to the effective mass in an FRW cosmology at finite temperature 
||10| , p!I[| , akin to the hard thermal loop resummation. Assuming that the heat bath arises 
from particles and radiation at high temperature and that at temperatures larger than the 
critical the system is in local thermodynamic equilibrium (LTE) with = 0, Einstein's 
equation for the Hubble 'constant' becomes 



HHt) = 



;i.6) 



'Pi 

with g being a constant that depends on the particle content of the theory. At the time 
at which the phase transition takes place, that is when T{tpt) = Tc, the Universe is still 
dominated by radiation, since oc m^/A^ >> m^/A for weakly coupled theories. Thus we 
reach the first conclusion: that if a phase transition is driven by the cooling of the expand- 
ing Universe, such a phase transition will occur during the radiation dominated era. The 
inflationary period of exponential (De Sitter) expansion will occur when T^{t) << m^/4A 
when the vacuum energy dominates, which will happen after several e-folds of radiation 
dominated expansion, once the effective temperature has red-shifted to almost zero. 

The small inhomogeneities of the CMBR restrict the scalar self coupling in these models 
to be A < 10^^^ (this will also be understood later when we compute the amplitude 
of scalar perturbations), and assuming LTE may be unjustified in these weakly coupled 
theories. However, even if LTE holds for short wavelength modes of the field theory it is 
unlikely to hold for the long wavelength modes. The reason for this is the following: at 
the phase transition when T(tpt) ?s Tc ~ m/y/X the Hubble constant is if ~ m(m/[XMpi\). 
Assuming the De Sitter stage of exponential inflation to occur at a GUT scale, this implies 
that mjy^l^ 10^^ Gev with the result that at the time of the phase transition dnitpt) « 
m~^, i.e. the horizon size is much smaller than the Compton wavelength of the scalar particle. 
This prevents thermalization of long wavelength modes even when the short wavelength 
modes may be strongly coupled to the bath and reach LTE. Thus our second conclusion: 
the phase transition will be strongly supercooled, with the long wavelength modes falling 
quickly out of LTE, even when LTE prevailed for short wavelength modes. The dynamics 
must necessarily be studied away from quasi-equilibrium and any approach based on effective 
potentials will miss important physics. 

Thus the phase transition can be described as a sudden quench from the high temperature 
phase at T > Tc into the low temperature phase with T << Tc on time scales much shorter 



than the microscopic time scales for thermalization and relaxation |T^. If the initial state 
was disordered in the sense that the expectation value of the order parameter is zero, then 
the phase transition occurs with the order parameter sitting at the top of the potential. 
Because of the symmetry (p — > —0 a state with = = will maintain vanishing value of 
the order parameter under the dynamics, and the rolling of the field down its potential hill 
necessary to end the inflationary stage must be understood as a consequence of the quantum 
fluctuations. 
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We have argued previously |T3|-[T5|] that under these circumstances long wavelength fluc- 
tuations will grow almost exponentially in a manner very similar to spinodal decomposition 
and phase separation in condensed matter systems. The quantum fluctuations as measured 
by the equal time two-point correlation function of the fleld < t) > must grow [p!^ -[T5 



in such a way that the mean square root fluctuation of the fleld 5(f) = y< ^^{x,t) > will 
eventually sample the minima of the potential and reach an equilibrium state. That is to 
say that at long times, when a quasi-equilibrium state has been achieved, it must be that 
>~ m^/A. Therefore the quantum fluctuations must become non-perturbatively 
large. Hence we reach a third conclusion: if the expectation value of the scalar fleld is at 
the top of the potential hill with vanishingly small velocity during the inflationary epoch, 
quantum fluctuations that are responsible for the process of phase separation will grow 
non-perturbatively large. Therefore the dynamics must necessarily be studied within a non- 
perturbative framework. 

Although new inflationary scenarios were previously studied |l8|j9|jT^, our work is rather 
different in that it addresses the description of the dynamics including self-consistently and 
non-perturbatively the non-equilibrium growth of quantum fluctuations and their effect on 
the dynamics of the scale factor. 



C. The puzzling questions: 

Having reached this conclusion on the non-perturbative growth of the quantum fluc- 
tuations one is presented with unsettling puzzles: typically the scalar fleld is written as 
t) = (j){t) + 6(f){x,t) where 0(t) is assumed to be the zero mode or expectation value 
of the scalar fleld and 6(j){x, t) the small quantum fluctuations that are ultimately respon- 
sible for metric perturbations. Inflation terminates when rolls down the potential hill 
and reaches the minimum, oscillating about it and eventually decaying into lighter particles 
leading to the reheating stage. 

However, the scenario that is envisaged here is that of a quenched or supercooled phase 
transition in which = throughout the evolution and the fluctuations grow to be very 
large and to eventually sample the minima of the potential. 

Therefore one is led to the questions: a) what is truly rolling down?, b) how does 
inflation end?, c) how sensitive is the dynamics to the initial value of the expectation value 
of the scalar field?, d) if the fluctuations grow to become non-perturbatively large will not 
they provide a large contribution to the energy momentum tensor and modify the FRW 
dynamics?, e) can one make sense of small fluctuations to calculate density perturbations?. 

Even when the initial value of 7^ 0, an inflationary stage requires that it be sufficiently 
small for the energy density to be dominated by the vacuum term. Under these circumstances 
the quantum fiuctuations will nevertheless grow, and the dynamics must address not only 
the rolling of the expectation value, but also the growth of spinodal fluctuations. 

The goal of the present work is to answer all of these questions. We will see below that 
it is precisely the large fluctuations that provide the sensible answers to all of the questions 
above, thus reconciling the naive picture of a scalar fleld rolling down the hill. Furthermore, 
a detailed analysis of the dynamics will allow us to quantify the answer to these questions 
and, in particular, will allow a profound interpretation of the effective zero mode, its initial 
conditions, and the small fluctuations that will provide the small density perturbations. 
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II. THE MODEL AND EQUATIONS OF MOTION 



Having recognized the non-perturbative dynamics of the long wavelength fluctuations 
we need to study the dynamics within a non-perturbative framework. We require that such 
framework be: i) renormalizable, ii) covariant energy conserving, iii) numerically imple- 
mentable. There are very few schemes that fulfill all of these criteria: the Hartree and the 
large approximation |£,13,n]. Whereas the Hartree approximation is basically a Gaussian 



variational approximation [|1^,0] that in general cannot be consistently improved upon, the 
large approximation can be consistently implemented beyond leading order [p!8|,p!9| and 



in our case it has the added bonus of providing many light fields (associated with Gold- 
stone modes) that will permit the study the effects of other fields which are lighter than 
the infiaton on the dynamics. Thus we will study the inflationary dynamics of a quenched 
phase transition within the framework of the large limit of a scalar theory in the vector 
representation of 0{N). 

We assume that the universe is spatially fiat with a metric given by eq. ( |1.1| ). The 
matter action and Lagrangian density are given by 



Sn 



d'^xa^{t) 



(2.1) 



vm = — 

^ ' 8N 



2NM^ 
A 



(2.2) 



7^ = 6 



'd{t) d^{ty 



(2.3) 



where we have included the coupling of $(x) to the scalar curvature TZit) since it will arise 
as a consequence of renormalization [^. 

The gravitational sector includes the usual Einstein term in addition to a higher order 
curvature term and a cosmological constant term which are necessary to renormalize the 
theory. The Lagrangian density for the gravitational sector is therefore: 



IGttG 2 ^ ' 



with K being the cosmological constant. In principle, we also need to include the terms 
afSfiu as they are also terms of fourth order in derivatives of the metric 
(fourth adiabatic order), but the variations resulting from these terms turn out not to be 
independent of that of TZ^ in the fiat FRW cosmology we are considering. 

The variation of the action ( p.l|) with respect to the metric g^i, gives us Einstein's equa- 
tion 



G 



8nG 



+ aH, 



Kg, 



-T, 



(2.4) 



where G^^, is the Einstein tensor given by the variation of y/—gTZ, Hnu is the higher order 



curvature term given by the variation of J—glZ 



and T^y is the contribution from the 
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matter Lagrangian. With the metric ( [L.l|) , the various components of the curvature tensors 
in terms of the scale factor are: 



GO = -3(a/a)2, (2.5) 

i//^ = -6 (^7^ + 3^v}j . (2.8) 

Eventually, when we have fully renormalized the theory, we will set = and keep as our 
only contribution to a piece related to the matter fields which we shall incorporate into 
T 



A. The Large N Approximation 



To obtain the proper large N limit, the vector field is written as 

$(x, t) = {a{x,t),7i{x,tj), 
with 7? an — 1-plet, and we write 

a{x,t) = VN<P{t) + x{x,t) ; {a{x,t)) = y/N<p{t) ; {x{x,t)) 



0. 



(2.9) 



To implement the large limit in a consistent manner, one may introduce an auxiliary field 
as in |jT9|. However, the leading order contribution can be obtained equivalently by invoking 
the factorization n,|20ll: 



X'^ 6(x^)x^ + constant , 
(vr • nf 2(7r')7f2 - (Tf^)^ + 0{1/N) 



^2 2 



{^')x' + n'{x'), 
{^')x. 



(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 



To obtain a large N limit, we define [n,pO| 



N-l 



7f(f,t) =tlj{x,t) {1,1, ■■■,!), 
where the large N limit is implemented by the requirement that 

{r)^Oil), (x^)^O(l), 0^O(1). 



(2.15) 



(2.16) 



The leading contribution is obtained by neglecting the 0{1/N) terms in the formal limit. 
The resulting Lagrangian density is quadratic, with linear terms in x and vr. The equations of 
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motion are obtained by imposing the tadpole conditions < xi^y t) >= and < n{x, t) >= 
which in this case are tantamount to requiring that the hnear terms in x ci-nd vf in the 
Lagrangian density vanish. Since the action is quadratic, the quantum fields can be expanded 
in terms of creation and annihilation operators and mode functions that obey the Heisenberg 
equations of motion 



7f(x, t) 



(Pk 

(2^ 



-ik-x 



The tadpole condition leads to the following equations of motion |14,20]: 



+ ?,H{t)^{t) + M\t)<p{t) = 0, 



with the mode functions 



fk{t) = 0, 



where 



M\t) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



In this leading order in 1/A^ the theory becomes Gaussian, but with the self-consistency 
condition 



d'k |/,(t)p 
(27r)3 2 



(2.21) 



The initial conditions on the modes fk{t) must now be determined. At this stage it 
proves illuminating to pass to conformal time variables in terms of the conformally rescaled 
fields 



see 



for a discussion) in which the mode functions obey an equation which is 
very similar to that of harmonic oscillators with time dependent frequencies in Minkowski 
space-time. It has been realized that different initial conditions on the mode functions lead 
to different renormalization counterterms in particular imposing initial conditions in 
comoving time leads to counterterms that depend on these initial conditions. Thus we chose 
to impose initial conditions in conformal time in terms of the conformally rescaled mode 
functions leading to the following initial conditions in comoving time: 



fkito) 



1 



fkito) 



d{to 



a{to 



fk{to), 



with 



W^ = e + M\to 
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(2.22) 



(2.23) 



At this point we recognize that when A^^(to) — TZ{to)/6 < the above initial condition 
must be modified to avoid imaginary frequencies, which are the signal of instabilities for 
long wavelength modes. Thus we define the initial frequencies that determine the initial 
conditions ( p.22| ) as 
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wl = k^ + 



M\to) - 



Wl = k^ + M\to)- 



6 
6 



for k < 



7^(^o 



for P > 



6 



(2.24) 
(2.25) 



In the large N limit we find the energy density and pressure density to be given by, 

r2\ 2 



e 1 -2, , A / 2M" 



+ - 



d^k 

2(2^ 



\m\'+u;i{tmt)[^ 



X 



(2.26) 



p + e 
N 



Sk 

2(27r) 



3a2(t) 



(2.27) 



It is straightforward to show that the bare energy is covariantly conserved by using the 
equations of motion for the zero mode and the mode functions. 

III. RENORMALIZATION 



Renormalization is a very subtle but important issue in gravitational backgrounds [£T . 
The fluctuation contribution {ip^{x,t)), the energy, and the pressure all need to be renor- 
malized. The renormalization aspects in curved space times have been discussed at length 
in the literature and has been extended to the large N self-consistent approximations 
for the non-equilibrium backreaction problem in |[1^J2^ . More recently a consistent and co- 
variant regularization scheme that can be implemented numerically has been provided (see 
J. Baacke's contribution to these proceedings and p^ , |23[| ). 

In terms of the effective mass term for the large N limit given by ( |2.2(JD and defining the 
quantity 



B{t) = a\t) (M\t)-TZ/Q 

A 



A 



B 



(3.1) 
(3.2) 



where the subscript B stands for bare quantities, we find the following large k behavior for 
the case of an arbitrary scale factor a{t) (with a(0) = 1): 



1 



+ 



+ 



ka'^{t) 2k^a^{t) 
1 



-B{t)] 



d r 



^^^^2 \B{t)[W{t)] + ait)- [a{t)B{t)\ \ + 0{l/e) 



(3.3) 



k 1 

+ 



H\t) + -(M\t)-TZ/Q 



a\t) ka'^{t) 

+ {-B{tf - a{tfi3{t) + ?>a{t)d{t)m - Ad\t)B{t)] + 0{llk% (3.4) 
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Although the divergences can be dealt with by dimensional regularization, this procedure 
is not well suited to numerical analysis (see however ref. [^). We will make our subtrac- 
tions using an ultraviolet cutoff constant in physical coordinates. This guarantees that the 
counterterms will be time independent. The renormalization then proceeds much in the 
same manner as in reference |Tl|; the quadratic divergences renormalize the mass and the 
logarithmic terms renormalize coupling constant and the coupling to the Ricci scalar. The 
renormalization conditions on the mass, coupling to the Ricci scalar and coupling constant 
are obtained from the requirement that the frequencies that appear in the mode equations 
are finite [O, i.e: 

- ml + ^snt) + ^<p\t) + ^{^\t))s = -ml + e^7^(^) + ^^^(t) + ^(^^(t))^ (3.5) 



2 ^ ' ' 2 

Finally we arrive at the following set of renormalizations 



2 
Of: 



StiGr 



StiG 



B 



2k«- 



1\ 



as 



Kb- 



m 



R 



m 



in - 

A^ 
167r2 
A2 



67 167r2 
2 ln(A//€) 



2 k/?- 



m 



ln(A/K) 
167r2 



R' 



167r2 
A2 



R 



Xr 

mt)) 



IB 



167r2 



167r2 



+ 



ml ln(A/K) 



167r2 



2 ln(A//s:) 



167r2 



ln(A/K) 
167r2 



Ar — A, 



ln(A//t) 
' 167r2 



2(27r) 



A;a2(t) 2W(t) 



7W2(t) -7^/6 



(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 
(3.11) 

(3.12) 



Here, k is the renormalization point. 



the renormalizations found using dimensional regularization [22,23 



As expected, the logarithmic terms are consistent with 

Again, we set aR = 



and choose the renormalized cosmological constant such that the vacuum energy is zero in 
the true vacuum. We emphasize that while the regulator we have chosen does not respect 
the covariance of the theory, the renormalized energy momentum tensor defined in this way 
nevertheless retains the property of covariant conservation in the limit when the cutoff is 
taken to infinity. In what follows, we drop the subscripts from the renormalized parameters. 

The logarithmic subtractions can be neglected because of the coupling A < 10"^^ Using 
the Planck scale as the cutoff and the infiaton mass as a renormalization point, these 
terms are of order X\n[Mpi/mR] < 10"^'', for m > 10^ GeV . An equivalent statement is 
that for these values of the coupling and infiaton masses, the Landau pole is well beyond the 
physical cutoff Mpi. Our relative error in the numerical analysis is of order 10~®, therefore our 
numerical study is insensitive to the logarithmic corrections. Though these corrections are 
fundamentally important, numerically they can be neglected. Therefore, in what follows, we 
will neglect logarithmic renormalization and subtract only quartic and quadratic divergences 
in the energy and pressure, and quadratic divergences in the fluctuation contribution. 
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IV. RENORMALIZED EQUATIONS OF MOTION FOR DYNAMICAL 

EVOLUTION 



It is convenient to introduce the following dimensionless quantities and definitions, 

Wu Xr 



H k 

T = niRt ; h = ; q = 

rriR rriR 



rriR 
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(4.1) 



2m% 



it) ; g^r) 



2m| 



fqir) = fk{t) 



(4.2) 



Choosing ^r = (minimal coupling) and the renormalization point k, = |m/j| and setting 
a(0) = 1, the equations of motion become: 



1^ + 3/.A_i +,^.) 



Vir) = 0, 



(4.3) 



f^(0) = [-hiO)-tco,]fM 

7^(o) 



Qitl'r 



q' + l-v\0) + '^-gm 



6m% 



for q^ > -1 +r]\0) 



for q^ < -1 +T]^{0) 



no) 

6m|. 



+ g^io), 



g^o). (4.4) 



The initial conditions for 7][t) will be specified later. An important point to notice is that 
the equation of motion for the q = mode coincides with that of the zero mode ( [4.31) . 
Furthermore, for //(r oo) 7^ 0, a stationary (equilibrium) solution of the eq. ( |4.3D is 



obtained when the sum rule 113114,20 



1 + 77^(00) + ^S(oo) = 



(4.5) 



is fulfilled. This sum rule is nothing but a consequence of Goldstone's theorem and is a 
result of the fact that the large approximation satisfies the Ward identities associated 
with the 0{N) symmetry, since the term —l + r]'^ + gH is seen to be the effective mass of 
the modes transverse to the symmetry breaking direction, i.e. the Goldstone modes in the 
broken symmetry phase. 

In terms of the zero mode //(r) and the quantum mode function given by eq.( [4.4| ) we find 
that the Friedmann equation for the dynamics of the scale factor in dimensionless variables 
is given by 



h\T) = 4hleR{T) 



hi 



^MIiXr 



(4.6) 
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and the renormalized energy and pressure are given by: 



2 
9_ 
2 



, , 2A^m|j f 1 .2 
(P + £)r = U"^ +9 I q dq 



l/.r-5^(?,r) + ^(|/,r-5(^)(g,r; 



(4.7) 
(4.8) 



where the subtractions S^^^ and S^'^^ are given by the right hand sides of eqns. (|3.4|) and 
respectively. 

We want to emphasize the following aspects of the set of equations ([4.3| , |4.4]j4r^ ): they 
determine the full dynamics of the matter plus classical gravity including backreaction effects 
both on the metric as well as in the dynamics of the fields, they are fully renormalized, 
maintain covariant conservation and in principle they can be consistently improved in the 
expansion. 

In order to provide the full solution we now must provide the values of ri{0), r]{0), and 
Hq. Assuming that the inflationary epoch is associated with a phase transition at the GUT 
scale, this requires that Nrnj^/Xji ^ (10^^ Gev )'^ and assuming the bound on the scalar 
self-coupling Xji ^ 10^^^ — 10^^^ (this will be seen later to be a compatible requirement), we 
find that Hq ^ N^^^ which we will take to be reasonably given by /lo ~ 1 — 10 (for example 
in popular GUT's ^ 20 depending on particular representations). 

We will begin by studying the case of most interest from the point of view of describing 
the phase transition: 77(0) = and ?7(0) = 0, which are the initial conditions that led to our 
puzzling questions. With these initial conditions, the evolution equation for the zero mode 
eq. (|4.3| ) determines that 77 (r) = by symmetry. 



A. Early time dynamics: 



Before engaging in a full numerical study, it proves illuminating to obtain an estimate 
of the relevant time scales and an intuitive idea of the main features of the dynamics. 
Because the coupling is so weak and after renormalization the contribution from the quantum 
fluctuations to the equations of motion is finite, we can obtain an estimate of the early time 
dynamics by neglecting the backreaction terms in the equations for the mode functions ( |4.4| ) 
and the Hubble constant ([4.6| ). Setting rj = and (7S f» in eq. ( [4.4|) and also setting to 
zero the terms proportional to g in eq.( |4.7| ), the evolution equations for the mode functions 
are those for an inverted oscillator in De Sitter space-time, which have been studied by Guth 



and Pi 



One obtains the approximate solution 



h{t) ^ ho, 




ho 



4 ' hr 



\ho 



(4.9) 
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where J±u{z) are Bessel functions, and Aq and Bg are determined by the initial conditions 
on the mode functions. 

When the physical wavevectors cross the horizon, i.e. when qe'^'^'^/hQ ^ 1 we find that 
the mode functions factorize: 

/.W-^^ (^ye(^-/^)'^-. (4.10) 

This result reveals a very important feature: u > 3/2; because of the negative squared mass 
term in the matter Lagrangian leading to symmetry breaking, we see that all of the mode 
functions grow exponentially after horizon crossing (for positive squared mass they would 
decrease exponentially after horizon crossing). This exponential growth is a consequence 
of the spinodal instabilities which is modified in De Sitter space-time but is a hallmark of 
the process of phase separation that occurs to complete the phase transition. We note, in 
addition that the time dependence is exactly given by that of the q = mode, i.e. the 
zero mode, which is a consequence of the redshifting of the wavevectors and the fact that 
after horizon crossing the contribution of the term 5'^/a^(r) in the equations of motion 
become negligible. Then we clearly see that the quantum fluctuations grow exponentially 
and they will begin to be of the order of the tree level terms in the equations of motion when 
5(E(r) ^ 1. At large times E(r) ^ J^'^{ho){ho/27iye^'^''-^^''°'^, with J='{ho) a finite constant 
that depends on the initial conditions and is found numerically to be of 0{1). 

In terms of the initial dimensionfull variables, this condition translates to < ip'^{x, t) >i?~ 
mj^/Xji, i.e. the quantum fluctuations sample the minima of the (renormalized) tree level 
potential. We find that the time at which the contribution of the quantum fluctuations 



becomes of the same order as the tree level terms is estimated to be 14 



At this time, the contribution of the quantum fluctuations makes the back reaction very im- 
portant and, as will be seen numerically, this translates into the fact that also determines 
the end of the De Sitter era and the end of inflation. The total number of e-folds during the 
stage of exponential expansion of the scale factor (constant Hq) is given by 

- In ( ttII^^I = Ihl In ( -1^1 + 0{l/h,) (4.12) 



(2z/-3) \\hlJ^^{ho) J 2" \\hlJ^^{ho 



For large we see that the number of efolds scales as as well as with the logarithm 



of the inverse coupling. These results ( |4.10| , ^?Tl| , [4.12| ) will be confirmed numerically below 



and will be of paramount importance for the interpretation of the main consequences of the 
dynamical evolution. 



B. r/(0) / 0: classical or quantum behavior? 

Above we have analyzed the situation when r]{0) = (or in dimensionfull variables 
0(0) = 0). The typical analysis of inflaton dynamics in the literature involves the classical 
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evolution of with an initial condition in which 0(0) is very close to zero (i.e. the top 
of the potential hill) in the 'slow-roll' regime, for which <^ 3-ff0. Thus, it is important to 
quantify the initial conditions on for which the dynamics will be determined by the classical 
evolution of and those for which the quantum fluctuations dominate the dynamics. We can 
provide a criterion to separate classical from quantum dynamics by analyzing the relevant 
time scales, estimated by neglecting non-linearities and backreaction effects. Considering 
the linear evolution of the zero mode in terms of dimensionless variables, and considering 
//(O) 7^ {ri{0) 7^ simply corresponds to a shift in origin of time), we find 



7]{t) ^ r]{0)e^''--^^''°\ (4.13) 

The non-linearities will become important and eventually terminate inflation when 77 (r) ^ 1. 
This corresponds to a time 



If Tc is much smaller than the spinodal time given by eq.( [4.11j ) then the c/asszca/ evolution 
of the zero mode will dominate the dynamics and the quantum fluctuations will not become 
very large, although they will still undergo spinodal growth. On the other hand, if Tc ^ 
the quantum fluctuations will grow to be very large well before the zero mode reaches the 
non- linear regime. In this case the dynamics will be determined completely by the quantum 
fluctuations. Then the criterion for the classical or quantum dynamics is given by 

r/(0) ^ VXho =^ classical dynamics 

ri{0) ^ \^ho =^ quantum dynamics (4-15) 

or in terms of dimensionfuU variables 0(0) ^ Hq leads to classical dynamics and 0(0) ^ Hq 
leads to quantum dynamics. 

However, even when the classical evolution of the zero mode dominates the dynamics, 
the quantum fluctuations grow exponentially after horizon crossing unless the value of is 
very close to the minimum of the tree level potential. In the large approximation the 
spinodal line, that is the values of for which there are spinodal instabilities, reaches all the 
way to the minimum of the tree level potential as can be seen from the equations of motion 
for the mode functions. Therefore even in the classical case one must understand how to 
deal with quantum fluctuations that grow and become large after horizon crossing. 



C. Numerics 



The time evolution is carried out by means of a fourth order Runge-Kutta routine with 
adaptive stepsizing while the momentum integrals are carried out using an 11-point Newton- 
Cotes integrator. The relative errors in both the differential equation and the integration are 
of order 10~*. We find that the energy is covariantly conserved throughout the evolution to 
better than a part in a thousand. Figures (1-3) show 5'S(r) vs r, /i(r) vs. r and ln(|/g(r)p) 
vs r for several values of q with larger q's corresponding to successively lower curves. Figures 
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(4,5) show p{T)/e{T) and the horizon size h~^{T) for A = 10~^^ ; 77(0) = ; rj{0) = and we 
have chosen the representative value Hq = 1.5. 

Figures 1 and 2 show clearly that when the contribution of the quantum fluctuations 
gT^ij) becomes of order 1 inflation basically ends, and the time scale for gT, to reach 0{1) 
is very well described by the estimate ( |4.11| ). From figure 1 we see that this happens for 
r 75, leading to a number of e-folds A^^e ^ 110 which is correctly estimated by (|4.11| , [4. 121) . 

Figure 3 shows clearly the factorization of the modes after they cross the horizon given 
by eq. ( |4.10|) . The slopes of all the curves after they become straight lines in figure 3 is 
given exactly by (2i/ — 3), whereas the intercept depends on the initial condition on the mode 
function and the larger the value of q the smaller the intercept because the amplitude of the 
mode function is smaller initially. Although the intercept depends on the initial conditions 
on the long-wavelength modes, the slope is independent of the value of q and is the same as 
what would be obtained in the linear approximation for the square of the zero mode at times 
long enough that the decaying solution can be neglected but short enough that the effect 
of the non-linearities is very small. Notice from the figure that when inflation ends and the 
non-linearities become important all of the modes basically saturate. This is also what one 
would expect from the solution of the zero mode: exponential growth in early-intermediate 
times (neglecting the decaying solution), with a growth exponent given by (z/ — 3/2) and 
an asymptotic behavior of small oscillations around the equilibrium position, which for the 
zero mode is r/ = 1, but for the 7^ modes depends on the initial conditions. All of 
the mode functions have this behavior once they cross the horizon. We have also studied 
the phases of the mode functions and we found that they 'freeze' after horizon crossing. 
This is natural since both the real and imaginary parts of fq obey the same equation but 
with different boundary conditions. After the physical wavelength crosses the horizon, the 
dynamics is insensitive to the value of q for real and imaginary parts and the phases become 
independent of time. Again, this is a consequence of factorization. 

The growth of the quantum fluctuations is sufficient to end inflation at a time given by Tg 
in eq. (|4.11| ). Furthermore figure 4 shows that during the inflationary epoch p{T)/e{T) ^ — 1 



and the end of inflation is rather sharp at Tg with p{T)/e{T) oscillating between ±1 with 
zero average over the cycles, resulting in matter domination. Figure 5 shows this feature 
very clearly; h{T) is constant during the De Sitter epoch and becomes matter dominated 
after the end of inflation with /i~^(r) ^ 3(r — Tf)/2, with Tf the time at the end of inflation. 
There are small oscillations around this value because both p(r) and e(r) oscillate. These 
oscillations are a result of small oscillations of the mode functions after they saturate, and 
are also a feature of the solution for a zero mode. 



D. Zero Mode Assembly 

This remarkable feature of factorization of the mode functions after horizon crossing can 
be elegantly summarized as 

fk{t)\kMt)«H=9iq,h)foiT), (4.16) 

with /o(t) a real function that obeys the zero mode equation ([4.3|), kph(t) = ke^^^ the 
physical momentum, and g{q, h) a complex constant. Since the factor g{q, h) depends solely 
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on the initial conditions on the mode functions, it turns out that for two mode functions 
corresponding to momenta ki, ^2 that have crossed the horizon at times ti > t2, the ratio of 
the two mode functions at time ts > t > ti > t2 is fkA^) I fkiif) oc e^''"!)''^^!"'^^) > x_ Then 
if we consider the contribution of these modes to the renormalized quantum fluctuations 
a long time after the beginning of inflation (so as to neglect the decaying solutions), we 
find that (7S(r) ^ (jg(2v-3)/ir _|_ gj^^all , where 'small' stands for the contribution of mode 
functions associated with momenta that have not yet crossed the horizon at time r, which 
give a perturbatively small (of order A) contribution. Then it is clear that after several 
e-folds from the beginning of inflation, we can define an 'effective zero mode' as 

Tjlffir) = gT,{T), or in dimensionfuU variables, 4>eff(t) = {4''^ , t)) ji ^ . (4-17) 

Although this identification seems natural, we emphasize that it is by no means a trivial or 
ad- hoc statement. There are several important features that allow an unambiguous identi- 
fication: i) [{ip'^ {x , t)) a] is a fully renormalized operator and hence finite, ii) because of the 
factorization of the superhorizon modes that enter in the evaluation of [{ip'^{x,t)) n], (f)eff{t) 
( [4.17] ) obeys the equation of motion for the zero mode, iii) this identification is valid sev- 
eral e-folds after the beginning of inflation, after the transient decaying solutions have died 
away and the integral in {ip'^{x, t)) is dominated by the modes with wavevector k that have 
crossed the horizon at t{k) << t. Numerically we see that this identification holds through- 
out the dynamics but for a very few e-folds at the beginning of inflation. This factorization 
determines at once the initial conditions of the effective zero mode that can be extracted 
numerically: after the first few efolds and long before the end of inflation we find 

</>e//W = 0e//(O)e(^-i)^* ; 0,^^(0) ^ ^^(/7/m), (4.18) 

where we parametrized 0e//(O) = ^JF(iJ/m) to make contact with the literature. We find 
numerically that J-'{H/m) ^ 0{1) for a large range of 0.1 < H/m < 50 and depends on the 
initial conditions of the long wavelength modes. 

The factorization of the superhorizon modes (that dominate the integral) implies that 

9jQ'dq\f^{r)\^C',\f',{r)\, (4.19) 
9j<l'dq\f^ir)\^C',\f',{r)\, (4.20) 

We also find numerically that even when ri{0) ^ this factorization phenomenon is very 
robust, and after a few e-folds from the beginning of inflation the dynamics is completely 
determined by the effective zero mode 



^e//(r) = V^2(r) + (?S(r). (4.22) 

We have checked numerically that the dynamics of the scale factor and equation of state 
obtained from the full quantum problem is exactly equivalent to that obtained from the 
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classical problem in terms of rjeffir). We have also checked numerically that the estimate 
for the classical to quantum crossover given by eq. (H:.15|) is quantitatively correct. Thus in 
the classical case in which 77(0) ^ we find that rjeffir) = r]{T) whereas in the opposite, 

quantum case Tjeffij) = ^gT^^r). We have run the numerical evolution of the scale factor 
with only the classical equation for r^g// with the proper initial conditions and found that 
it coincides within our numerical error with the evolution obtained by the full system of 
equations in either the classical or quantum case. 

This remarkable feature of the zero mode assembly of long-wavelength spinodally un- 
stable modes is a consequence of the presence of the horizon. It also explains why despite 
the fact that asymptotically when the fluctuations sample the broken symmetry state, the 
equation of state is that of matter. Since the excitations in the broken symmetry state are 
massless Goldstone bosons one would expect radiation domination. However, the 'assembly' 
phenomenon, i.e. the redshifting of the wave vectors, makes these modes behave exactly like 
zero momentum modes that give an equation of state of matter domination (upon averaging 
over the small oscillations around the minimum). 



V. MAKING SENSE OF 'SMALL FLUCTUATIONS': 

Having recognized the effective classical variable that can be interpreted as the compo- 
nent of the field that drives the FRW background and rolls down the classical potential hill, 
we want to recognize unambiguously the small fluctuations. We have argued above that 
after horizon crossing, all of the mode functions evolve proportionally to the zero mode, 
and the question arises: which modes are assembled into the effective zero mode and which 
modes are treated as perturbations? In principle every k ^ mode provides some spatial 
inhomogeneity, and assembling these into an effective homogeneous zero mode seems in prin- 
ciple to do away with the very inhomogeneities that one wants to study. However, scales 
of cosmological importance today have first crossed the horizon during the last 60 or so 
e-folds of inflation. Recently Grishchuk |P3 has argued that the sensitivity of the measure- 



ments of AT/T probe inhomogeneities on scales ~ 500 times the size of the present horizon. 
Therefore scales that are larger than these and that have first crossed the horizon much 
earlier than the last 60 e-folds of inflation are unobservable today and can be treated as an 
effective homogeneous component, whereas the scales that can be probed experimentally via 
the CMB inhomogeneities today must be treated separately as part of the inhomogeneous 
perturbations of the CMB. 

Thus a consistent description of the dynamics in terms of an effective zero mode plus 
'small' quantum fluctuations can be given provided: a) the total number of e-folds N^. ^ 60, 
b) all the modes that have crossed the horizon before the last 60-65 e-folds are assembled 

into an effective classical zero mode via 4>eff(t) = [0o(^) + (^^^(2^5 ^))-r] ^ 5 c) the modes that 
cross the horizon during the last 60-65 e-folds are accounted as 'small' perturbations. The 
reason for the requirement a) is that in the separation t) = (f)eff{t) + 6(f){x, t) one requires 
that 5(f)[x,t)/(j)eff{t) ^ 1. As argued above after the modes cross the horizon, the ratio of 
amplitudes of the mode functions remains constant and given by e('^-i)^^ with AN bemg 
the number of e-folds between the crossing of the smaller k and the crossing of the larger 
k. Then for t) to be much smaller than the effective zero mode, it must be that the 
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Fourier components of 5(j) correspond to very large fc's at the beginning of inflation, so that 
the effective zero mode can grow for a long time before the components of 50 begin to grow 
under the spinodal instabilities. In fact requirement a) is not very severe; in the figures (1-5) 
we have taken = 1.5 which is a very moderate value and yet for A = 10~^^ the inflationary 
stage lasts for over 100 e-folds, and as argued above, the larger ho for fixed A, the longer 
is the inflationary stage. Therefore under this set of conditions, the classical dynamics of 
the effective zero mode 0e//(^) drives the FRW background, whereas the inhomogeneous 
fluctuations 6(j){x, t) , which are made up of Fourier components with wavelengths that are 
much smaller than the horizon at the beginning of inflation and that cross the horizon during 
the last 60 e-folds, provide the inhomogeneities that seed density perturbations. 



VI. SCALAR METRIC PERTURBATIONS: 



Having identified the effective zero mode and the 'small perturbations', we are now in 
position to provide an estimate for the amplitude and spectrum of scalar metric pertur- 
bations. We use the clear formulation by Mukhanov, Feldman and Brandenberger [25| in 



terms of gauge invariant variables. In particular we focus on the dynamics of the Bardeen 
potential p6|, which in longitudinal gauge is identified with the Newtonian potential. The 
equation of motion for the Fourier components (in terms of comoving wavevectors) for this 
variable in terms of the effective zero mode is lESl 



H{t) 



+ 2 H{t) - H{t) 



$fc = 0. 



(6.1) 



We are interested in determining the dynamics of for those wavevectors that cross 
the horizon during the last 60 e-folds before the end of inflation. During the inflationary 
stage the numerical analysis suggests that to a very good approximation 



H{t) ^ Ho ; <Peff{t) = 0e//(O)e(^-^)^°*, 
where Ho is the value of the Hubble constant during inflation, leading to 



Ho 



Ho 



(3 = u-l 



(6.2) 



(6.3) 



The coefficients 0^,6^ are determined by the initial conditions. 

Since we are interested in the wavevectors that cross the horizon during the last 60 
e-folds, the consistency for the zero mode assembly and the interpretation of 'small pertur- 
bations' requires that there must be many e-folds before the last 60. We are then considering 
wavevectors that were deep inside the horizon at the onset of inflation. Mukhanov et. al. 
25| show that ^k{t) is related to the canonical 'velocity field' that determines scalar per- 



turbations of the metric and which is quantized with Bunch-Davies initial conditions for the 
large fc-mode functions. The relation between and v and the initial conditions on v lead 
at once to a determination of the coefficients and for k » Ho 



25 



ak = -- 



6TT 



(0). 



TT 



■f 



2Hok 



0. 



(6.4) 
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Thus we find that the amphtude of scalar metric perturbations after horizon crossing is 
given by 



Mt)\ = kl\Mt)\^l 



(0) 



2Ho 
k 



(6.5) 



The power spectrum per logarithmic k interval is given by |5fc(t)p. The time dependence of 
\Sk\ displays the unstable growth associated with the spinodal instabilities of super-horizon 
modes and is a hallmark of the phase transition. This time dependence can be also under- 
stood from the constraint equation that relates the Bardeen potential to the gauge invariant 



field fluctuations , which in longitudinal gauge are identified with S(j){x, t) . To obtain 
the amplitude and spectrum of density perturbations at second horizon crossing we use the 



conservation law associated with the gauge invariant variable 



6 



1 



0, 



(6.6) 



which is valid after horizon crossing of the mode with wavevector k. Using this conservation 
law, and the relation that during the inflationary stage 1+p/e = 8n4>ljj-/2MpiHQ <^ 1, and 
assuming matter domination at second horizon crossing and $fc(^/) = p5[' ^^"^ 



Mtf)\ 



12r(z/)V? /2Ho 



5(z/ - l)J^{Ho/m) V k 



(6.7) 



where J-'{Ho/m) determines the initial amplitude of the effective zero mode ( [4.18| ). We can 
now read the power spectrum per logarithmic k interval 



leading to the index for scalar density perturbations 

1 - 2(u 



2' 



(6.8) 



(6.9) 



We remark that we have not included the small corrections to the dynamics of the effective 
zero mode and the scale factor arising from the non-linearities. These are expected to lead 
to perturbatively small 0{\) corrections to the index (|0|). The spectrum given by ( |6.7| ) 
is similar to that obtained in references PJTO[] although the amplitude differs from that 
obtained there. We emphasize an important feature of the spectrum: it has more power 
at long wavelengths because z/ — 3/2 > 0. This is recognized to be a consequence of the 
spinodal instabilities that result in the growth of long wavelength modes and therefore in 
more power for these modes. This seems to be a robust prediction of new inflationary 
scenarios in which the potential has negative second derivative in the region of field space 
that produces inflation. 

It is at this stage that we recognize the consistency of our approach for separating the 
composite effective zero mode from the small fluctuations. We have argued above that many 
more than 60 e-folds are required for consistency, and that the 'small fluctuations' correspond 
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to those modes that cross the horizon during the last 60 e-folds of the inflationary stage. For 
these modes H^/k = e~^''**(^) where t*{k) is the time of horizon crossing of the mode with 
wavevector k since the beginning of inflation. The scale that corresponds to the Hubble 
radius today Aq = 27i/ko is the first to cross during the last 60 or so e-folds before the 
end of inflation. Smaller scales today will correspond to k > ko at the onset of inflation 
since they will cross the first horizon later and therefore will reenter earlier. The bound 
on \6ko\ oc AT/T < 10~^ on these scales provides a lower bound on the number of e-folds 
required for these type of models to be consistent: 

12 ln(z/-i) 
iV,>60 + -— ^-^-3^, (6.10) 

where we have written the total number of e-folds as Ne = iJo^*(^o) + 60. This in turn can 
be translated into an upper bound on the coupling constant using the estimate given by 
eq.(m- 



VII. CONTACT WITH THE RECONSTRUCTION PROGRAM: 

The program of reconstruction of the inflationary potential seeks to establish a rela- 
tionship between features of the inflationary scalar potential and the spectrum of scalar and 
tensor perturbations. This program, in combination with measurements of scalar and tensor 
components either from refined measurements of temperature inhomogeneities of the CMB 
or through galaxy correlation functions will then offer a glimpse of the possible realization of 
the inflation 1^^^. Such a reconstruction program is based on the slow roll approximation 
and the spectral index of scalar and tensor perturbations are obtained in a perturbative 
expansion in the slow roll parameters P7|J2^ 



V{<P) = -4t ■ (7-2) 
H(p 

We can make contact with the reconstruction program by identifying above with our 0e// 
after the first few e-folds of inflation needed to assemble the effective zero mode from the 
quantum fluctuations. We have numerically established that for the weak scalar coupling 
required for the consistency of these models, the cosmologically interesting scales cross the 
horizon during the epoch in which H ^ Hq ; (peff ~ {i^ ^3/2) Hq (p^ff V ~ m%/ X ^ 4'lff- 
In this case we find 

V{ci>eff) = ^ - ^ ; <^eff) ^ 0{\) « r/(0e//). (7.3) 



With these identifications, and in the notation of p7| , |28| the reconstruction program 
predicts the index for scalar density perturbations Us given by 

n.-l = -2(z/-^) + 0(A), (7.4) 
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which coincides with the index for the spectrum given by the power spectrum per logarithmic 
interval with \6k\ given by eq. (|6.7|) . We must note however that our treatment did not 
assume slow roll for which (z/ — |) -C 1. Our self-consistent, non-perturbative study of the 
dynamics plus the underlying requirements for the identification of a composite operator 
acting as an effective zero mode, validates the reconstruction program in weakly coupled 
new inflationary models. 

VIII. CONCLUSIONS: 

We have studied the non-equilibrium dynamics of a new inflation scenario in a self- 
consistent, non-perturbative framework based on a large N expansion, including the dy- 
namics of the scale factor and backreaction of quantum fluctuations. Quantum fluctuations 
associated with superhorizon modes grow exponentially as a result of the spinodal insta- 
bilities and contribute to the energy momentum tensor in such a way as to end inflation 
consistently. 

Analytical and numerical estimates have been provided that establish the regime of 
validity of the classical approach. We find that these superhorizon modes re-assemble into 
an effective zero mode and unambiguously identify the composite operator that can be used 
as an effective expectation value of the inflaton field whose classical dynamics drives the 
evolution of the scale factor. This identification also provides the initial condition for this 
effective zero mode. 

If the model allows many more than 60 e-folds we provide a criterion that allows the 
identification of the small perturbations that give rise to scalar metric (curvature) perturba- 
tions. We then use this criterion combined with the gauge invariant approach to obtain the 
spectrum for scalar perturbations. We find the index to be less than one, providing more 
power at long wavelength as a result of the spinodal instabilities. We argue that this 'red' 
spectrum is a robust feature of potentials that lead to spinodal instabilities in the region in 
field space associated with inflation. Finally we made contact with the reconstruction pro- 
gram and validated the results for these type of models based on the slow-roll assumption, 
despite the fact that our study does not involve such an approximation. A more detailed 
version of this article with a discussion of the issue of decoherence is forthcoming p9| . 
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FIG. 3. ln(|/g(r)|2) vs. r, for r?(0) = 0,?7(0) = 0, A = IQ-^^/io = 1.5 for q = 0.5,2.0,5.0,8.0,10 
smaller q corresponds to larger values 
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